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S ! Abstract 

The J/ifjDD form factor is evaluated in a QCD sum rule calculation for both 
D and J/ip off-shell mesons. We study the double Borel sum rule for the 
CN ■ three point function of two pseudoscalar and one vector meson current. We 

find that the momentum dependence of the form factors is different if the D 
or the J/ifi meson is off-shell, but they lead to the same coupling constant in 
the J/^DD vertex. 
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The coupling constant of the strong interactions at the most fundamental level, i.e., 
of quarks and gluons, has been studied both theoretically in lattice QCD calculations and 
Ph. experimentally in e + e~ colliders. Our knowledge on a s has significantly increased during 
the last decade |]. Unfortunately, these advances did not bring a better understanding of 
the seemingly more accessible quantities: the coupling constants bewteen hadrons. Indeed, 
on this next level of the strong interactions there has been little progress in what regards 
the determination of the "hadron strong charges". Even in the most studied case, the 
NNir coupling, there is still some controversy. In the strange sector, the couplings between 
nucleons, hyperons and strange mesons have been strongly constrained by the use of SU(3) 
symmetry and intense phenomenological analyses of low energy hadronic reactions 0. In 
the charm sector, neither it is reasonable to use SU(4) symmetry since it is badly broken, nor 
experimental information is available. Of course, there are, here and there, exceptions as, 
for example, the D*D7r coupling, which was recently determined by the CLEO collaboration 
from D* decays ||. 

In a curious development of the recent history of hadron physics, the couplings involv- 
ing charm mesons became a very important building block in the construction of effective 
lagrangian theories designed to describe low energy reactions such as, for example: 

7T + J/tfj -> D* + D (1) 

Actually, apart from the couplings, more detailed knowledge is required, namely, we need to 
know the form factors at these effective vertices, as for example gD*Dn(Q 2 )- These effective 



1 



theories may be used to compute cross sections between light and charmed mesons, which, 
in turn may be used to understand charm hadron production (and suppression) in the 
nucleus-nucleus collisions performed at RHIC [f4|-|i"3| . 



The work of calculating these form factors is rewarded with an interesting side product: 
information about the size of charm mesons. Of course, the size of a hadron depends on how 
we "look" at it. The most extensively studied particle is the nucleon, which has been probed 
mainly by photons. In lower energy experiments, where also the four momentum transfer (q 2 ) 
is low, it was possible to determine the electromagnetic form factor (and the charge radius) 
of the nucleon. In higher energies experiments and very large values of Q 2 (Q 2 = —q 2 ) a very 
different picture of the nucleon emerged, in which it is made of pointlike particles, the quarks. 
From these observations one may conclude that, when probing the nucleon, nearly on-shell 
photons (q 2 ~ 0) recognize sizes whereas highly off-shell photons (q 2 << 0) do not. This 
statement is supported by the phenomenologically very successful vector meson dominance 
hypothesis, according to which real photons are with a large probability converted to vector 
mesons (which are extended objects) and then interact with the nucleon. 

A couple of years ago we started our program of computing the above mentioned quan- 



tities in the framework of QCD sum rules (QCDSR). In |L4| we calculated the coupling 



constant in the vertex NDA C with the D meson off-shell. In |TJJ we did a similar calculation 



for strange coupling constants. In [16| we extended the calculation performed in [II and 



computed the Q 2 dependent form factors of the NDA C and ND*A C vertices. We have also 
studied the D*Dir, B*Bir 0,0] and DDp |19[ vertices. One of the conclusions of these 



works is that when the off-shell particle in the vertex is heavy, the form factor tends to be 
harder as a function of Q 2 , which means larger cut-off parameters and smaller associated 
sizes. 

In the present work we will further investigate form factors involving heavy mesons in 
order to extend our previous conclusions. We also want to better estimate the uncertainties 
in the procedure of determining coupling constants with our techniques. For these purposes 
we consider the vertex J/ipDD and compute form factors and coupling constants for the 
cases where the D is off-shell, the J/ip is off-shell and then compare the results. 

Following the QCDSR formalism described in our previous works we write the three-point 
function associated with a J/ tpDD vertex, which is given by 



TfW) = J d 4 x d^y(0\T{j D (x)jUy)jtm\0)^- x e- i ^-y , (2) 
for a D off-shell meson, and by 

T^W) = / d'xdS^\T{ 3 D{x)jl{y)Um^)e ip, - x e-^-^ , (3) 

for a J/ij; off-shell meson, where jr> = idjsc, y D = icr/^d and = c'j^c are the interpolating 
fields for the outgoing D + , incoming D + and J/ip respectively with d and c being the down, 
and charm quark fields respectively. 

The general expression for the vertex function, or correlator, in Eqs. (|^) and (^) has two 
independent structures. We can write T M in terms of the invariant amplitudes associated 
with these two structures: 
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r M (p,p') = ri(p 2 ,p ,2 ,g 2 )p M + r 2 (p 2 ,^,g 2 w 



l2 



(4) 



As usual we shall write the correlators (|J) and (||) in terms of hadron and quark-gluon 
degrees of freedom. These two representations, also called phenomenological and theoretical 
side, are then indentified one with the other yelding a sum rule. On the theoretical side we 
make an operator product expansion (OPE), keeping the first terms, which are represented 
in Fig. la and lb for the correlators (Q) and (|3]) respectively. The diagram of Fig. la 
and the first of Fig. lb may be evaluated giving the perturbative contributions. Applying 
Cutkosky's rule to these contributions we can write a double dispersion relation for each one 
of the invariant amplitudes Y\ ert (i = 1,2), over the virtualities p 2 and p' 2 holding Q 2 = —q 2 
( with q = p' — p) fixed: 
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where pf ert (s,u J Q 2 ) equals the double discontinuity of the amplitude Ti(p 2 ,p' 2 , Q 2 ) on the 
cuts s min < s < oo, m 2 < u < oo, and where s min = 4m 2 in the case of the D off-shell, where 
the dispersion relation is written in terms of the two D mesons' momenta, and s m i n = m 2 in 
the case of J / ip off-shell, the dispersion relation being now written in terms of the J / ip and 
the D meson momenta. In terms of the s, u and t invariants, the double discontinuities are: 



pert(D) 
Pi 
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+ (s-t-w + 2m 2 )(l 
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for a D off-shell meson, with the integration limits for u: 

1 



u 



max, mm 



2m 2 



st + mj(s + 2t) ± J(s - 4m 2 c )s(m 2 c - tf 



and 



(8) 



pPert(J/tp) 



A3/2 



[s — t — u)(m^ — su) 



(9) 



pert(J/ip) 
P2 



A 3 / 2 



(s + t — u)(rric — su) 



for a J/ip off-shell meson, with the corresponding integration limits for u: 

1 



max,mtn 



2ml 



-st + m 2 (2s + t)±(s- mi)J-t(Am 2 - t) 



(10) 



(11) 



with A = X(s,t,u) = s 2 + t 2 + u 2 — 2st — 2su — 2tu and t = —Q 2 . For this last limit we 
remark that the condition u > t — m 2 must be satisfied and that, as it can be seen from the 
square root, in the timelike region, t = q 2 > 0, we must have t > 4m 2 . 
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Since we are dealing with heavy quarks, we expect the perturbative contribution to be 
dominant on the OPE side. However, for the J /if) off-shell it turns out that the contribution 
of the quark condensate, shown in the second diagram of Fig. lb, is also important and must 
be included: 

T<^W) = tf-^-rrfl ^ + ti = r<??> (Pm + P^) (12) 

The phenomenological side of the vertex function is obtained by considering the contri- 
bution of the J/ if) and one D meson, or the two D mesons states to the matrix element in 
Eqs. (|2|) and ([3]) respectively. In doing so, we introduce the decay constants fo and fj/^, 
which are defined by the matrix elements 

(0\3d\D) = ^ , (13) 
m c 

and 

Wb» = mj/^e; , (14) 
where e M is the polarization of the vector meson. We also make use of the matrix element: 
(J/if;(q)D(p')\D(p)) = - gj/ ^ DD (Q 2 )( Pfl+ p'^ . (15) 
Altogether these steps lead to: 

T^(D) M = _?&D f {Q2) 1 1 x 

p m c p* — rrij,^ p — m D 

m 2 D + m 2 /tf! + Q 2 \ 
—2p H ^ Pn + higher resonances , (lb) 

m J/i> j 



Tf en ^\p,p') = -^Mkgj^Q 2 ) 1 1 U + p\ + higher resonances , (17) 

f m z it n z — rn.%, TV — mi, V ^ ' 



l pfp 

m 2 ' p* — m D p'~ — m D 

where 



unU?) " >,hl " ^ • (is) 

and 



2 r J D ) 
D 

m c Q 2 + tti 2 d 



>/d 9ddj/^(Q 2 



gj/AQ 2 ) = wjtf/ /° 2 ^ m ? , (19) 



where m D , mj^ are the masses of the mesons D and J /if) respectively, and g^pj/^iQ 2 ) is 
the form factor at the DDJ/ip vertex when the meson M is off-shell. The contribution of 
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higher resonances in Eqs. ([16]) and (|PT|) will be taken into account in the standard form of 
continuum contribution from the thresholds Sq and Uq J20 . 



For consistency we use in our analysis the QCDSR expressions for the decay constants 
appearing in Eqs. (|l6f) and fllTD up to dimension four: 



fl 
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it 
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f-(qq)e<- m B- 

D 



■ml) /M„ 



(20) 



i no (r + 2m 2 c )Jr-4ml , , _, /M 2 



where u = (m D + A u ) 2 GeV 2 , r = (rrij/^ + A r ) 2 GeV 2 and M m 2 is the Borel mass used in 
the two-point functions given above. 

We have omitted the numerically insignificant contribution of the gluon condensate. 

Eq. (|T7| ) shows that when the J/ip is off-shell, the sum rules in both structures give the 
same results. On the other hand, according to Eq. fllfi|), when the D is off-shell the two 
structures give different results. In this last case we will concentrate in the p'^ structure, 
which we found to be the more stable one. The sum rules in this structure read: 

T£ ' >(p ,p ,q ) = T y 2 ' >(p ,p ,q ) (22) 

for a D off-shell and 

■y^pheni J/ip) , 2 /2 2\ -r^pertU / ij}) / 2 /2 2\ , -p<W>/ 2 /2 2\ /oo\ 

r 2 lv \p , P , q ) =r 2 {p ,p ; ? ) + r \p ,p , <? ) (23) 

for & J/ip off-shell. 

Inserting in these equations the corresponding expressions for the perturbative, quark 
condensate and phenomenological terms and performing a double Borel transformation |2^ 



in both variables P 2 = —p 2 — > M 2 and P' 2 = — p' 2 — >• M' 2 we obtain the final expressions for 
the sum rules. In order to allow for different values of M 2 and M' 2 we take them proportional 
to the respective meson masses, which leads us to study the sum rules as a function of M 2 
at a fixed ratio: 

Ml = !% (24) 



M' z m D 



for flU) and M 2 = M' 2 for 

In refs. [2~I,Z2]] it was found that relating the Borel parameters in the two- (M m 2 ) and 
three-point functions (M 2 ) as 

2M 2 m = M 2 , (25) 

is a crucial ingredient for the incorporation of the HQET symmetries, and leads to a con- 
siderable reduction of the sensitivity to input parameters, such as continuum thresholds sq 
and Mo, and to radiative corrections. Therefore, in this work we will use Eq. (^) to relate 
the Borel masses. 
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The parameter values used in all calculations are m q = (m u + m^)/2 = 7 MeV, m c = 
1.3 GeV, m D = 1.87 GeV, m J/t = 3.1 GeV, (qq) = -(0.23) 3 GeV 3 . We parametrize the 
continuum thresholds as 

s = (m M + A s ) 2 , (26) 

where uim = mclmj/^) for the case that the J/ip(D) meson is off-shell, and 

u = (m D + A u ) 2 . (27) 

Using A s = A u = A r ~ 0.5 GeV for the continuum thresholds and fixing Q 2 we found 
good stability of the sum rule for 9^ooj/^{Q 2 ) as a function of M 2 in the interval 4 < M 2 < 

16 GeV 2 and also for g^oJ/^Q 2 ) in the interval 2 < M 2 < 10 GeV 2 . 

Fixing M 2 = 11 GeV 2 we calculate the momentum dependence of the form factor 
^DDj/iiiQ 2 ) m ^ ne interval —0.5 < Q 2 < 5.0 GeV 2 , where we expect the sum rules to 
be valid (since in this case the cut in the t channel starts at t ~ m 2 and thus the Euclidian 
region stretches up to that threshold). Our numerical calculations can be well reproduced 
by the gaussian parametrization: 

&»(Q 2 ) = 16-4 e-' 12 ^' (28) 

We stress here that it was not possible to fit our results with a monopole form! 
As in ref. 0, we define the coupling constant as the value of the form factor at Q 2 = 
-rnL, where mu is the mass of the off-shell meson. In the case of an off-shell D, this leads 



to 



9DDJ/* = 8 - 05 ( 29 ) 

Along the same lines we choose M 2 = 8 GeV 2 and calculate the momentum depen- 
dence of the form factor g^Dj/^Q 2 )- The numerical results can be fitted by the monopole 
parametrization: 

9ddjUQ 2 ) = q 2 + 143 . 18 ( 3 °) 
which, at the J/ip pole leads to: 

9%% = 7-98 (31) 

All our numerical results and parametrizations are shown in Fig. 2. The circles and 
squares correspond to the numerical results for the D and J/ip off-shell respectively. These 
points are fitted by the dashed (Eq. (|28|) ) and solid lines (Eq. (|30|) ) respectively and 
extrapolated to the meson poles which are represented by the triangles, which give the 
numbers quoted in (|29|) and (j3l"l). 

A closer look into Fig. 2 reveals that, whereas the circles are well fitted by the dashed 
line, this is not the case of the squares, which are not very accurately described by the solid 
line. Due to this reason and also because the J/ij) pole is farther away, the estimate fl3"T| ) 
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is much less reliable than (|29|). In fact, to reduce the uncertainty in the extrapolation of 
the J/ip off-shell form factor, we have to impose the condition that both coupling constants 
dr^Dj/ip an d 9ddj/iP mus t coincide. We therefore use the former as a guide in the fitting 
procedure (and subsequent extrapolation) leading to the latter. This condition also imposes 
a severe constraint in the minimum Q 2 used in the calculation. We have used the interval 
5.5 < Q 2 < 9.5 GeV 2 . Including smaller values of Q 2 in the calculation would increase the 
curvature of the solid line, which, when extrapolated to lower Q 2 values, would give coupling 
constants very different from g^j/^- In other words, it would become impossible to obtain 
the two triangles at the same height. 

The star in Fig. 2 indicates the J/ip off-shell form factor taken at Q 2 = 0. As it can be 
immediately obtained from ([ft]), g^DJ/^ = 7.47 at Q 2 = 0. This coincides with the estimate 
made in Ref. [|J using the vector meson dominance model. 

The information presented in Fig. 2 may now be crossed with those obtained in [|19|| , 
about the size of the D meson. In Fig. 3 we compare 9ddj/^{Q 2 ), given by ([ftj) with 
9dd p (Q 2 )i calculated in |19| and parametrized by: 



9^ ] Dp {Q 2 ) = 2-53 e-& (32) 

The DDJ/ip form factor is represented by a solid line whereas the DDp one is shown with 
a dashed line. 

On a qualitative level, the comparison between solid and dashed lines both in Fig. 2 and 
Fig. 3 shows that the form factor is harder if the off-shell meson is heavier. In particular, 
remembering that there is an overlap between the photon and the vector mesons p and J/ip, 
we check in Fig. 3 an empirical formula regarding the resolving power of the photon, used 
some time ago by experimentalists. HERA data on electron-proton reactions could be well 



understood introducing a "transverse radius of the photon", parametrized as [23]: 



r 7 ~ 1/ VQ 2 + m 2 (33) 

where m is the mass of the vector meson considered. This empirical formula tells us that for 
Q 2 — > oo the photon is pointlike and "resolves" the nucleon target, i.e., identifies its pointlike 
constituents and does not "see" the size of the nucleon. Moreover this formula indicates that 
for Q 2 ~ and for light mesons (like the p°) the photon has appreciable transverse radius 
and therefore also identifies the global nucleon extension. Finally, in the above formula we 
may have a heavy vector meson (J/ip or T) which will, either real or virtual, resolve the 
nucleon into pointlike constituents. This conjecture can be applied to a D target probed by 
a J/^and seems to be supported by Fig. 3 when we consider the vicinity of Q 2 ~ 0. There, 
we see (again) that the slope of the dashed line is stronger than that of the solid line. This 
means that also around the "real photon" region, the J/ip resolves smaller scales then the 
P- 

In Fig. 4 we compare g'oDj/^iQ 2 )-' gi ven by (f2"£|) (solid line), with g^)D P {Q 2 ) (dashed line), 



calculated in [19] and parametrized by: 
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We can clearly observe that the D + p — > D transition has a much harder form factor than 
the D + J ftp -> D one. 

In the limiting case rj/^ « r D << r p these results can be understood in a simple 
picture: when the D hits the smaller J/ip it can "see" a size, it can measure it! On the 
other hand, when the D hits the much larger p meson, it does not see any size, in the same 
way as large Q 2 photons (in DIS measurements) do not see any size in the proton. Rather, 
they will interact with pointlike partons! In real life the mentioned radii are not so different 
from each other and therefore the differences in the curves in Fig. 4 are not so pronounced. 

Our calculations contain uncertainties in the QCD parameters (masses and vacuum con- 
densates), in the OPE (because we neglect higher order operators), in the model of the 
continuum (the uncertainty in the values of so and uq) and the systematic error in the 
extrapolation procedure to obtain the couplings. These errors are present in most of the 
QCDSR calculations and are to a certain extent unavoidable. Due to them most of QCDSR 
results are plagued by a 20% error. Therefore our numbers for the couplings have this same 
uncertainty. On the other hand, the curves shown in the figures are so dramatically different 
that they confirm beyond any doubt our previous suspiction [p^|IT| , |T9"f regarding transition 
vertices and, most of all, as one can see clearly in Fig. 3, they show in the case of "charge 
form factors" (where the same meson goes in and out the vertex), that the D meson "seen" 
by a p is much larger than when it is probed by the J/ip. 

Our program will continue and studies of other vertices are in progress. 
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Q 2 (GeV 2 ) 

FIG. 2. Momentum dependence of the DDJip form factor. Circles and squares represent our 
numerical calculations for the D and J/ip off-shell respectively. The dashed and solid lines give 
the parametrization of the QCDSR results through Eq. (f2§|) for the circles and Eq. (|30|) for the 
squares. The triangles give the form factors at the poles of the particles (which we indentify with 
the coupling constant). The star shows the form factor at Q 2 = 0. 
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Q 2 (GeV 2 * 

FIG. 3. Momentum dependence of the DDJif; (solid line) and DDp (dashed line) form factors. 
In both cases the vector mesons are off-shell. 




Q 2 (GeV 2 ) 



FIG. 4. Momentum dependence of the DDJt/j (solid line) and DDp (dashed line) form factors. 
In both cases the D meson is off-shell. 
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